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ABSTRACT 

We  consider  a  controlled  system 

X  =  f(X/U) 

exhibiting  a  self  excited  periodic  solution  x(t)  for  u(t)  =  0  and  consider 
the  question  of  modification  of  this  orbit  to  a  controlled  periodic  solution 
x(t)  corresponding  to  a  control  u(t)  chosen  so  as  to  minimize  a  cost 
functional  of  the  form 

A 

“  w(x(t),  u(t) )dt  , 

T 

A 

where  T  is  the  period  of  the  controlled  periodic  solution*  Some  relevant 
applications  are  cited. 
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SI6MIFICAMCE  AND  BXPIAHATION 


Great  pains  have  been  taken  to  make  Section  1  of  this  paper  an 
explanation  of  the  background  and  significance  of  the  work.  This  is 
particularly  true  of  the  first  three  pages  of  the  section.  This  material 
requires  minimal  technical  background  and  is  recommended  to  the  reader  wishing 
to  appreciate  the  basic  ideas  of  the  paper  without  getting  too  involved  in  the 
technical  material. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


OPTIMAL  ORBITAL  REGULATION  IN  DYNAMICAL  SYSTEMS 
SUBJECT  TO  HOPE  BIFURCATICm* 

** 

David  L*  Russell 

1.  Introduction.  Many  systems  Involving  self-excited,  nonlinear  oscillations 
can  be  usefully  modelled  by  a  system  of  the  form 

X  -  A{e)x  +  B(6)  [J  +  D(e)u  ,  (1.1) 

y 

X  e  E*',  USE™  , 

y  +  ■Y(®#y)y  +  lc{6)y  -  c(fl)  x  +  o(6)  u  ,  (1.2) 

y  scalar. 

The  first  system  (1.1)  In  which  A(d),  B(6),  d(8)  are  matrices  of  dimension 
n^n,  nx2,  n’^m,  respectively,  represents  a  linear  oscillator  of  some 
sort  with  a  degree  of  Internal  damping,  l.e.,  the  eigenvalues  of  A(8)  have 
negative  real  parts  for  all  values  of  the  parameter  B  under  consideration. 
The  scalar  second  order  equation  (1.2)  in  y  represents  a  nonlinear 
oscillator.  In  It  c(B),  e(6)  are  vectors  of  dimension  n,  m,  respectively, 
and  *  denotes  transpose.  In  most  cases  the  uncoupled  equation 

y  +  Y(9,y)y  +  k(e)y  -  O  (1.3) 

has  y  >  y  ■  0  as  an  asymptotically  stable  critical  point  for  B  In  some 

range,  say  0  <  B  <  8^,  that  critical  point  becoming  unstable  for  ®  ^  ®q» 

•  • 

Typically  the  term  Y(B,y)y  talses  the  form 
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■V(8,y)y  -  y^y  -  Y^(fl)y  +  Y(0»y)y  , 

where  >  0,  9y^(8)  >0,  6  f*  0.  The  higher  order  teme  represented  by 

A  *  *  * 

Y(Bry)y  provide  strong  restoring  forces  when  y  is  large.  The  familiar 
consequence,  amply  documented  in  the  literature  (see,  e.g.,  [A],  [B] ,  [C] )  is 
that  for  ^  (1.3)  has  a  stable  periodic  solution  expanding  rapidly  away 

from  the  origin  (0,0)  as  0  Increases  beyond  8^. 

The  uncontrolled  coupled  system  ((1.1),  (1.2)  with  u  =  0)  typically 
involves  a  matrix  B(6)  which  is  rather  small,  so  that  the  matrix  describing 
the  linearization  at  the  origin, 

A(0)  B(e) 

is  nearly  lower  block  triangular  and  is  a  stability  matrix  for  8  <  8^, 
becoming  unstable  for  8  >  8^ ,  where  8^,  the  bifurcation  point,  is  near  the 

8^  value  already  discussed  for  (1.3).  The  system  (1.1),  (1.2)  then  likewise 
exhibits  a  stable  periodic  solution  for  8  >  8^,  lying  close  to  the  y,  y 
plane  for  8  near  8^, 

One  of  the  simplest  examples  is  the  Hartlen-Curry  model  ( [D] )  for  aero- 
dynamically  Induced  flutter  of  a  long  elastic  rod  of  circular  cross  section  as 
shown  schematically  in  Figure  1.1.  A  stream  of  fluid,  e.g.,  air,  flo«rs  past 
the  rod  with  velocity  V.  For  a  certain  "flutter  speed",  Vg,  the  fluid  flow 
is  steady  for  V  <  Vq.  For  V  >  Vg  the  fluid  flow  is  no  longer  steadyi 
alternating  vortices  form  in  the  wake  behind  the  rod,  producing  a  near- 
periodic  force  of  alternating  direction  with  a  frequency,  dependent  on  V, 
which  is  known  as  the  Strouhal  frequency.  The  Hartlen-Curry  model  uses 
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equations  of  the  form  (cf.  [B] ) 


•  2 

X  +  ftt  +  Px  “  aV  y 

y  +  (o  -  aV)y  +  ~  y^  +  fiv^y  ■  bx 

to  describe  the  self-excited  oscillations  of  the  coupled  system  for  V  >  Vq 
(Vq  is  slightly  greater  than  0/a),  Typically  3,  o,  a,  b  are  rather 
small  quantities. 


Figure  1.1  Oscillations  of  a  rod  with 
circular  cross  section. 

A  more  complicated  example,  for  which  we  will  not  write  out  the  equations 
in  detail,  is  of  considerable  interest  in  flight  dyneunics.  It  is  rather 
similar  to  the  example  just  presented  for  the  rod  with  circular  cross  section, 
except  that  we  now  envision  a  wing  cross-section  immersed  in  the  flowing  fluid 
as  shown  in  Figure  1.2.  The  displacements  shown  in  the  figure  are,  of  course, 
exaggerated  by  comparison  with  actual  operating  levels  (one  hopesl). 


Figure  1*2 

Here  z  represents  the  vertical  displacement  of  the  wing  tip,  t  is  the 
torsional  displacement  of  the  wing  at  the  tip,  and,  again,  y  is  the  vortex 
strength  measured  at  the  trailing  edge  of  the  wing*  One  obtains  here  a  six 
dimensional  system 

••  •  •  • 

Z  ■  f (z,Z,T,T,y,y,u,V) 

**  •  •  • 

T  -  g(z,Z,T,T,y,y,u,V) 

**  •  •  • 

y  -  h{z,Z,T,T,y,y,u,V) 

exhibiting  self  excited  oscillations  much  as  in  the  Hartlen-Curry  case  with 
added  features,  such  as  divergence  of  the  periodic  solution  to  infinity  at 
certain  parameter  values  for  V* 

The  initial  oscillations  arising  in  systems  of  this  sort  are,  at  least 
insofar  as  the  x  component  (cf*  (1.1),  (1*2))  is  concerned,  of  rather  small 
amplitude*  The  exception  which  causes  greatest  concern  occurs  when  the 
nonlinear  oscillator  equation  (1*3)  has  a  periodic  solution  whose  frequency  is 
close  to  one  of  the  natural  frequencies  of  vibration  of  the  linear  elastic 
system  modelled  by  x  ■  A(9}x*  When  this  situation  obtains,  quite  dramatic 
Increases  in  amplitude  may  result  in  (1*1),  (1*2)*  Figure  1*3  shows, 
qualitatively,  the  sort  of  results  that  one  obtains  with  a  Hartlen-Curry  type 
model* 
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Figure  1.3 

Among  other  things#  control  may  be  used  in  (1.1)#  (1.2)  to  increase  the 
range  of  parameter  values  9  for  which  the  origin  is  asymptotically  stable# 
l.e.#  to  increase  6^  and/or  to  modify  the  nature  of  the  oscillations  Which 
do  occur  for  6  >  9^  -  by  suppressing  the  amplitude  of  those  oscillations#  for 
example • 

The  traditional  mode  of  control  is  linear  feedback#  the  general  form  of 
which#  for  (1.1)#  (1.2)  would  l>e 

u  ■  Kx  +  l(,)  #  (1.5) 

y 

possibly  with  gains  scheduled  so  that  K  «  K(9),  L  *  L(9).  xt  will  ]9e 
recognised#  however#  that  direct  measurement  of  y#  y  will  often#  perhaps 
usually#  hm  impractical.  Certainly  the  measurement  of  the  vortex  strength  at 
the  trailing  edge  of  a  wing  would  require  a  sensory  tour  de  force.  So#  in 
practice#  (1.5)  must  ordinarily  Ise  modified  to 
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u  *  Kx  .  (1*6) 

When  one  investigates  how  (1.6)  affects  (1.4)/  one  sees  that  the  result  is  to 
replace  that  matrix  by 


A(8)  +  0(6)  K 


c(8)*  +  e(e)*K 


B(0) 


[-V»l  -To 


0 


(1.7) 


One  may  verify  that  If  T  satisfies  the  quadratic  matrix  equation 

AT  -  tF  -  TCT  +  B  -  0  (1.8) 

then  (1.7)  Is  similar  to 


A  -  TC 
C 


CT 


°)  . 

*  r/ 


(1.9) 


For  ^  ^  ^  unstable  and  it  may  be  assumed  that  X  Is  selected  so  that 

A  remains  stable.  It  Is  then  classical  that 

AT  -  tF  +  B  -  0  ^  (1.10) 

has  a  unique  solution  T  and  an  easy  application  of  the  implicit  function 
theorem  shows  that  for  small  B  (1.8)  has  a  unique  solution  near  the  solution 
of  (1.10).  Moreover f  an  B  tends  to  sero,  T  tends  to  sero.  Wow  our 
capability  to  influence  the  eigenvalues  of  CT  t  F  lies  entirely  in  the 
term  CT  appearing  in  (1.9).  If  B  is  quite  small/  large  changes  in  K  are 
required  to  produce  modest  changes  in  T.  It  follows  that  when  B  is  small, 
as  is  conmonly  the  case,  it  will  be  difficult  to  materially  affect  the 
eigenvalues  which  are  responsible  for  the  bifurcation  j^anomena,  l.e.,  the 
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eigenvalues  of  CT  +  T.  This  means  that  restricted  feedback  of  the  form  (1*6) 
Is  not  likely  to  be  very  effective  in  extending  the  range  of  values  of  0  for 
which  the  origin  remains  asymptotically  stable* 

It  may  be  argued  that  feedback  of  the  form  (1.5)  may  be  very  nearly 

realized  through  the  use  of  a  state  estimator  ( [F] )  for  y«  y.  However,  it  is 

likely  that  the  estimation  will  be  very  difficult  in  practice  because,  when 
« 

B  is  near  0,  y,  y  are  nearly  unobservable  via  measurements  on  x  and, 
moreover,  the  oscillations  of  (1.2)  xvecessarlly  take  place  in  a  region  where 
the  nonlinear  terms  appearing  there  balance  the  linear  terms,  "energy  wise", 
and  hence  are  of  comparable  magnitude.  Then  one  Is  trying  to  construct  a 
state  estimator  for  nearly  unobservable,  nonlinear  phenomenal  a  rather  heroic, 
not  to  say  Slsyphlc,  task. 

It  follows  that  if  linear  feedback  is  to  be  used.  It  will  take  the  form 
(1.6)  and  the  main  benefit  will  be  realized  throu^  choosing  K  so  that  the 
term  B(y)  (,]  In  (1.1)  affects  x  minimally  In  some  appropriate  sense.  For 

y 

example,  the  eigenvalue  of  A(B}  4-  D(9)X  corresponding  to  the  fundamental 
mode  of  vibration  can  be  moved  further  to  the  left  in  the  complex  plane  or  Its 
imaginary  part  can  be  Increased  so  that  resonance  occurs  at  higher 
frequencies,  thus,  in  effect,  pushing  the  "spike"  of  Figure  1.3  to  the  right. 
The  main  point,  insofar  as  our  current  discussion  Is  concerned.  Is  that  one 
must.  In  effect,  concede  that  oscillation  is  going  to  take  place  and  then  take 
steps  to  modify  or  suppress  the  manifestations  of  such  oscillations  In  the 
system  represented  by  (1.1). 

Having  seen  that  we  may  as  well  admit  from  the  first  that  a  self-excited 
periodic  motion  will  be  present  in  the  operating  system,  we  may  as  well  cite 
this  as  an  excuse  and  an  opportunity  to  deal  directly  with  such  solutions  as 
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we  attempt  regulation  of  the  ayetem  (1.2)*  Our  purpose  In  this  paper 

is  to  discuss  the  optimal  regulation  of  periodic  solutions  of  systems  of  this 
typ«f  primarily  with  a  view  toward  amplitude  suppression. 

Zn  Sections  2,  3  we  study  a  general  system  (wherein  6  Is  suppressed 
until  needed) 

X  *  f(x«u)  (1.11) 

for  X  e  Of  em  open  set  In  1^/  and  u  6  E^.  Most  of  our  results  are 
actually  obtained  In  the  context  of  the  specialized  systeois  of  the  form 

X  g(x)  +  H(x)u  . 

Our  objective  Is  to  show  the  existence  of,  and  to  characterize/  In  terms  of 
necessary  conditions,  state  and  control  pairs  x,  u  such  that  x  Is  periodic 
with  period  T  and  such  that  x,  minimize/  relative  to  an  appropriate 
family  of  state  and  control  pairs  x.  u  (x  with  period  T).  a  cost 
functional  of  the  general  form 

J(X/U/T)  ■  ^  w(x(t)/  u(t))dt  .  (1.12) 

In  Section  2  we  begin  by  supposing  that  the  uncontrolled  system 

X  -  f(x/0) 

has  a  periodic  solution  with  period  Tq  and  then  obtain  some  local  existence 
results  ensuring  the  existence  of  an  optimal  control,  at  least  In  an 
appropriate  local  sense.  Then  In  Section  3  we  Indicate  the  form  of  the 
necessary  conditions,  some  computational  methods  and  describe  some 
computational  work  already  undertaken,  concluding  with  an  Indication  of  work 
In  progress  and  envisioned  for  the  future. 


2<  Local  Bxletence  of  a  Solution  for  the  Periodic  Optimal  Control  Problem 


In  our  study  of  existence  questions  we  will  confine  attention  to  systems 
In  which  the  control  appears  linearly < 

X  ••  f(x,u)  ■  g{x)  +  H(x)u  ,  (2.1) 

wherein  It  Is  assumed  that,  for  some  open  subset  0 C  , 

n  nm 

g  :  0  R  ,  H  t  0  +  r'““ 

are  defined  and  continuously  differentiable  throughout  0«  Here  denotes 

the  space  of  n  x  m  matrices  with  real  entries.  The  differentiability  of  f 
with  respect  to  both  x  and  u  for  x  0  O,  u  0  sP,  Is  then  clear. 

Our  basic  assumption  Is  that  the  "uncontrolled"  system 

X  -  f(x,0)  ■  g(x)  (2,2) 

has  a  solution  x(t)  which  Is  periodic  with  least  positive  period  Tq, 

Beyond  that  we  are  concerned  with  the  variational  system  based  on  x(t), 
which  Is 

5  -  G(x(t))5  +  H(x(t) )u  ,  (2,3) 

G(x(t))  =  (x(t))  (Jacobian  of  g  w,r,  to  x) , 

9x 

We  denote  by  4(t,s)  the  fundamental  matrix  solution  of 

*  G<x(t))»(t,s)  ,  (2,4) 

4(s,s)  ■  1  (n  xn  Identity  matrix)  (2,5) 

and  abbreviate  4(t,0)  by  9(t),  Both  (2,3)  and  (2,4)  may  be  taken  to  be 
defined  for  all  real  t  If  we  extend  x(t)  In  the  obvious  way  by 
periodicity.  The  resulting  systems  are  periodically  dependent  on  t  with 
period  Tg,  The  period  transition  map  associated  with  the  periodic  solution 
x(t)  is  the  map  defined  by  the  matrix 


Since  the  system  (2.2)  Is  autonomous,  we  may  stipulate  any  point  on  the 
solution  x(t)  as  the  Initial  point  Xq  «  x(0).  The  assumption  that  x(t) 
is  a  non-trivial  periodic  solution  enables  one  to  see  readily  that 

p^  -  x(0)  -  ^  0  .  (2.6) 

It  is  well  known  that  p^  is  an  eigenvector  of  corresponding  to  the 

eigenvalue  ^  ^(T^).  This  is  an  immediate  consequence  of  the  fact 

that 

5(t)  s  Mt)  5  ♦(t)p^ 

• 

is  a  periodic  solution  of  C  ■  G(x(t))C.  We  assume  that  the  remaining  n  -  1 

eigenvalues  of  ^CT^),  which  we  denote  by  different  from 

1.  An  important  special  case  arises  when 

|X^|  <  1,  i  -  2,3, ...,n  ,  (2.7) 

in  which  case  the  periodic  solution  x(t)  of  (2.2)  is  locally  asymptotically 

stable.  The  condition  (2.7)  is  not  needed  in  this  section  but  becomes  very 

important  in  numerical  considerations  to  be  introduced  later. 

We  denote  by  {p^}  the  one  dimensional  subspace  of  R**  spanned  by  p^ 

and  by  P  the  (n-1 )-dimensional  subspace  spanned  by  the  (possibly 

generalized)  eigenvectors  P2,  P3/..*,P,^  of  ^(T^)  associated  with  the 

eigenvalues  X  ,  X  , ...,X  .  if  p  is  a  real  n  (n-1)  matrix  of  rank 
2  3  n 

n-1  whose  columns  are  P2,  or  consist  of  real  and  imaginary 

parts  of  P2,  P3,..«/P„  in  the  case  of  complex  X^,  X^,...,X^,  then  each 
pep  has  the  unique  representation 

p  -  po,  o  e  r""^  . 

Let  g^  be  the  unique  vector  in  R**  such  that 
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1,  s,  u  -  1  , 


p  “  o»  pep  . 

* 

Then  ie  an  eigenvector  of  4(7^)  corresponding  to  the  eigenvalue  1> 

We  denote  by  {q^}  the  one  dimensional  subspace  of  spanned  by  q^  and 

by  Q  the  (n-1 )-dlmenalonal  subspace  of  spanned  by  the  eigenvectors 

* 

*^2'  ^3'***'*ln  <^BSOclated  with  the  eigenvalues  of  *(Tq)  (If 

we  were  to  allow  «(T)  to  be  complex,  would  be  replaced  here 

by  1.,  X,,  •  •  • ,  X  )  •  Letting  Q  be  the  n  x  (n-1)  matrix  whose  columns  are 
2  J  n 

?2'  *l3'***'?n  real  and  Imaginary  parts)  It  may  be  arranged  that 

Q*P  -  in-1  • 


Our  requirement 


X^  1,  1  •  2,3,»»»,n  , 

guarantees  that  with 

e(t)  -  Q*(*(t)  -  1  )P  ,  (2.8) 

n 

0(Tg)  Is  a  nonsingular  (n-1)  x  (n-1)  matrix.  Evidently  we  have  the 
decompositions 

{p^ }  «  P  -  r"  -  {q^ }  +  Q  . 

Having  taken  care  of  the  various  definitions  and  assumptions  needed,  we 
turn  now  to  consideration  of  periodic  solutions  of  the  controlled  system  (2.1) 
corresponding  to  controls 

u  e  b^IO.T^],  T,  >  Tjj  , 

or,  more  precisely,  to  restrictions  of  such  controls  as  explained  below.  The 
theorem  which  we  present  below  may  be  proved  much  more  simply  when  we  require 
that  u  8  C[0,T^1  by  Invoking  the  implicit  function  theorem  and, 
additionally,  a  local  uniqueness  result  is  obtained  in  that  case.  Because  of 
the  method  of  proof  employed,  and  In  the  Interests  of  more  colorful 
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math«matlcal  metaphor,  we  ehall  call  thla  the  "arrow  <md  target"  theorem*  It 
will  turn  out  that  the  proof  rests  on  the  Intermediate  value  theorem. 

Theorem  ,2.1.  For  arbitrarily  small  positive  x,  6,  x  satisfying 
X  <  T^  -  Tg,  there  exists  e  ■  e( x, 6)  >  0  such  that  whenever 

e 

lul  <6  (2.9) 

B  1 


the  system  (2.1)  has  at  least  one  solution  w(t)  satisfying 

lw(t)  •  x(t)  I  <  S,  0  <  t  <  T^  ,  (2.10) 

and  such  that 

w(T)  -  w(0)  (2.11) 

for  at  least  one  value  of  T  with 

It  -  TqI  <  X  .  (2.12) 

Once  w,  T  corresponding  to  u  have  been  identified,  w,  u  may  be 
regarded  as  a  periodic  state,  control  pair  by  first  restricting  w  and  u 
to  (0,T]  and  then  defining 

w(t  +  KT)  -  w(t),  u(t  +  kT)  -  u(t),  k  -  ±1,  ±2,*».  .  (2.13) 

Proof  of  Theorem  2. 1 .  For  simplicity  we  begin  with  the  equation 

w(t)  -  g(w(t))  +  H(r(t))u  (2.14) 

wherein  r(t)  is  a  given  continuous  n- vector  valued  function  defined  on 
[0,T]  and  satisfying 


lr(t)  -  x(t)  I  „  <  6,  t  e  tO,T,]  , 
_n  1 


x(t)  being  the  periodic  solution  of  (2.1)  corresponding  to  u(t)  =  0  which 
we  have  discussed  above.  This  Includes,  of  course,  the  special  case  wherein 
H(x(t))  =  H,  H  constant  n  x  m  matrix  . 
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Once  we  have  proved  the  existence  of  a  periodic  pair  w,  u#  for  this  special 

case  we  will  be  able  to  describe  rather  easily  the  modifications  required  to 

prove  the  corresponding  result  for  (2.1). 

2 

Let  u  e  L  [0,T  ]  satisfy  (2.9)  with  e  yet  to  be  determined,  and  let 
in  1 

w(t)  be  the  solution  of  (2.1)  corresponding  to  the  control  u  and  the 
initial  state 

w(0)  »  X-  +  p,  p  e  P,  Ipl  <  p  ,  (2.15) 

°  r" 

for  some  p  >  0  yet  to  be  determined.  It  is  an  easy  consequence  of  the 
Caratheodory  existence  and  regularity  theory  (see  [6],  e.g.)  that  we  can 
assume 

lw(t)  -  x(t)  I  <  6  (2.16) 

r" 

A 

for  any  6  >  0  by  choosing  e  and  p  in  (2.9),  (2.15)  to  be  sufficiently 

small.  Our  intent  is  to  keep  u  fixed  emd  show  that  we  cam  find  p  as  in 

(2.15)  and  T  satisfying  (2.12)  such  that  (2.11)  holds. 

Let  y( t)  be  the  solution  of 

y(t)  -  g(y(t))  +  H(r(t))u(t)  ,  (2.17) 

y(0)  -  Xq  .  (2.18) 

Clearly  y(t)  satisfies  (2.16),  i.e., 

Iy(t)  -  x(t)  I  „  <  ?  (2.19) 

R 

if  e  is  sufficiently  small.  Then 

w(t)  -  y(t)  +  z(t)  (2.20) 

where 

z(t)  -  g(y(t)  +  z(t))  -  g(y(t))  ,  (2.21) 

z(0)  -  p  .  (2.22) 

The  variational  equation  associated  with  the  solution  y(t)  is 


C(t)  -  GCy(t))C(t) 


(2.23) 


wher«  G(y(t))  «  ^  (y(t})  la  tha  Jacobian  of  g  avaluatod  along  tha 

A 

solution  y«  If  m  lat  ^(t)  ba  tha  natrlx  solution  of 

♦(t)  -  G(y(t))*(t),  t  e  tO.T^l  ,  (2.24) 

♦(0)  -  I  (2.25) 


and  set 

«(t)  -  Q*(^(t)  -  DP  ,  (2.26) 

the  usual  regularity  arguments  show  that  8(t}  converges  unlfomdy  to  0(t)/ 

given  by(2.^).  as  lu*  .  and  hence  sup  {ly(t)  -  x(t)l  })  tends 

I.^[0,T^1  te[.0,T^)  r" 

to  zero.  Since  ^(T^)  has  been  seen  to  be  nonsingular.  It  follows  that 

A  A 

6(T)  Is  Invertible  with  6(T)  satisfying  some  bound 

I9(T)"^I  <B  ,  (2.27) 

provided  (2.9)  and  (2.12)  hold  with  e  and  t  sufficiently  snail. 

The  "target"  of  the  "arrow  and  target"  theorem  Is  the  Image.  Z^.  of  the 
subspace  P  c  r"  under  the  map 

Z(p.T)  -  z(p.T)  -  p  . 

where  z(p.t)  «  z(t)  la  the  solution  of  (2.21)  corresponding  to  the  Initial 
state  (2.22).  l.e..  z(p.O)  «  z(0)  *  p.  We  can  obtain  a  parametrized 
representation  of  z^  by  writing 

z(p.T)  -  P  ■  o^p^  +  p,  p  e  P  .  (2.28) 

•  ^ 

*  x(0)  •  g(x(0))  as  noted  earlier  and  (since  p  6  P)  setting 

®2 

P  -  pa.  o  -  (•  )  e  r"‘^  .  (2.29) 

a 

n 


with  P  the  n  ^  (n-1)  matrix  described  earlier.  Then 

* 

a^  -  q^(z(p.T)  -  p)  . 
a  »  Q*(z(p.T)  -  p) 


(2.30) 

(2.31) 
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and  writing 


P  -  p$,  e  e  r"“^  , 

we  see  that  the  Jacobian  of  ^  with  respect  to  3  at  3  »  0  (i.e*,  p  -  0) 

A 

is  precisely  6(T)  as  given  by  (2.26)*  From  this  and  the  implicit  function 
theorem  it  follows  that  can  be  represented,  for  p  near  0  and  T  near 

Tg,  as  corresponding  to  the  set  of  points  (2.30),  (2.31)  for  which 

-</(o,T)  , 

with  of  class  near  o“0,  T“Tq.  In  fact  (2.30)  shows  that 


(o,Tjj)  “  gp  (z(p,Tq)  -  pnp,o  “Illt^ck^O 


q*(«(T^j)  -  I)P§(Tg)“^ 


(2.32) 


The  surface  Z,^  will  thus  be  nearly  tangent  to  the  stibspace  P  at  p  >  0 

for  small  lul,  i.e.,  for  y(t)  near  x(t),  and  T  near  Tg,  because 

* 

♦(Tg)P  ■  PA  for  an  appropriate  (n-1)  x  (n-1)A  and  q^P  *  0. 

The  "arrow"  of  our  theorem  is  the  vector  Xg  -  y(T).  The  "arrow"  hits 

the  "target"  just  in  case  there  is  a  vector  p  e  P  and  an  instant  T  such 

that 

-  y(T)  -  y(0)  -  y(T)  -  z(p,T)  -  p 

-  z(p,T)  -  z(p,0)  -  z(T)  -  z(0)  , 


for  then 

w{0)  •  y(0)  +  z(0)  -  y(T)  +  z(T)  ■  w(T)  .  (2.33) 

That  such  a  "hit"  takes  place  is  a  consequence  of  the  intermediate  value 
theorem  for  continuous  functions  on  R^.  The  situation  in  hand  is  represented 
graphically  in  Figure  2.1  where  *•■*■•*  is  the  trajectory  of  Xg  >  x(T) 
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tor  T  near  Tq,  >  »■  -»  la  the  trajectory  of  Xq  -  y(T),  the 

hyperplane  P  is  shown,  and  the  curves  -----  Indicate  the  outline  of  a 
cylindrical  region 

be  -  (Xq  -  x(T))l  <  d  ,  (2.34) 

for  d  <  5,  truncated  by  hyperplanes  parallel  to  P 

H*  «  {x  I  X  -  -Yp^  +  p,  p  e  P,  p^x(-  • 


Figure  2. 1 

The  surface  depends  on  u  via  the  dependence  of  z  on  y  In 

equation  (2.21).  Indicating  this  dependence  by  Z,p(u),  we  begin  by  taking 
d  small  enough  In  (2.34)  so  that  Z^^(O)  bisects  the  cylinder  Into  two 
regions.  This  must  be  the  case  for  d  small  because  the  tangent  vector  to 
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Xg  -  x(T)  at  T  «  Tg  la  p^,  i^ila  the  tangant  hyperplana  to  ^^(0)  at 
the  origin  is  P*  Next,  L  and  T  are  restricted  so  that 
p*(Xg  -  x(T))  >  2L,  Tg  -  2T  <  T  <  Tp  -  T  , 

P*(Xq  -  X{T))  <  -21,,  Tjj  +  T  <  T  <  Tg  +  2T  , 

Then  we  make  £,  and  hence  lul  ,  small  enough  so  that  Xg  -  y(T) 

m  1 

lies  in  the  cylinder  (2«34)  for  |T  -  Tg|  <  T  and 

pJ<Xg  -  y(T))  >  L,  Tg  -  2T  <  T  <  Tg  -  T  ,  (2.35) 

p*(Xg  -  y(T))  <  -L,  Tg  +  T  <  T  <  Tg  +  2T  ,  (2.36) 

both  of  %rtiich  are  possible  because  y(t)  converges  uniformly  to  x(t)  as 
lui  .  tends  to  sero.  Finally,  e  is  further  restricted,  if 

li^(0,Ti] 

necessary , 

so  that  z,n(u)  bisects  the  cylinder  (2.34)  when  lul  ,  <  t.  Because 

z(0,T)  3  0,  2(0, T)  =  0  and  the  surface  Z,p(u)  always  passes  throu^  the 
3  ip 

origin.  Then  —  (p,*!)  can  be  uniformly  bounded  for  p  small  and  T  near 
3a 

Tg,  using  (2.27),  (2.32),  and,  further  restricting  d  if  necessary,  the 
intersection  of  Z,p(u)  with  the  cylinder  (2.34)  does  not  meet  the  ends,  R^, 
of  the  cylinder  for  lul  _  <  6,  |t  -  T-|  <  T,  Put  another  way,  if 

P^(T)  denotes  the  connected  neighborhood  of  the  origin  in  P  which  maps  into 
the  cylinder  (2.34)  under  Z(*,T), 

|p*2(p,T)l  <  L  for  lul  2  <  e.  It  -  Tg |  <  T,  p  e  P^(T)  .  (2.37) 

In  I 

It  follows  that  we  may  unambiguously  denote  the  two  components  of  ths  cylinder 
(2.34)  cut  out  by  Z,p  by  C'*'(T),  C~(T),  according  to  whether  that  component 
contains  H***  or  H*,  respectively.  For  x  in  the  cylinder  (2.34)  and 
Tg  -  2T  <  T  <  Tg  ♦  2T  we  define 

D(x,T)  -  ±  dls(x,Z^)  , 
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th«  diatano*  balag  tha  dlatanca  within  tha  cyllndar  (2.34)  and  *  balng  usad 
In  C'*’(T),  -  In  C~(T).  Than  conaldar  tha  function  D(Xg  ~  y(T).T).  Worn  tha 
continuity  of  D(x.T)  (tha  proof  of  which  wa  omit  hara>  but  It  la  not  raally 
difficult  to  aatabllah)  It  followa  that  0(xq  -  y(T)fT)  la  a  contlnuoua 
function  of  T.  From  (2.35),  (2.36)  It  followa  that  D(Xq  -  y(T),T)  la 
poaltlva  for  -  2t  <  T  <  -  t,  nagatlva  for  +  t  <  T  <  +  2t.  From 

the  Intarmadlata  value  theorem  there  axlata  at  laaat  one  T  In 
Tq-t<T<Tq  +  t  auch  that  Dg(xQ  -  y(T),T)  ■  0,  from  which  It  followa 
that 

Xg  -  y(T)  e  Zj 

for  auch  T  and,  aa  Indicated  earlier,  thla  ahowa  that  there  axlata  w(t)  a 
aolutlon  of  (2.14)  auch  that 

w(0)  •  w(T) 

and  we  have  the  reault  for  the  ayatem  (2.14). 

We  Indicate  only  briefly  how  the  reault  la  extended  to  the  ayatem  (2.1), 

1  .e . , 

w  ■  g(w)  +  H(w)u  .  (2.1r) 

The  main  complication  la  that  when  we  let 

y  -  g(y)  +  H(y)u,  y(0)  -  ,  (2.38) 

then.  In  order  that  w  •  y  +  z  ahould  be  a  aolutlon  of  (2.1),  we  muat  have 

*  -  g(y+a)  -  g(y)  +  (H(y+a)  -  H(y))u  (2.39) 

ao  that  amall  aolutlona  s(t)  of  (2.21)  do  not  approximately  aatlafy  (2.23)* 
but  rather  the  equation 

5"  tG(y(t))  +  0^(t,y(t) )]  c  (2.40) 


where 


Is  an  n  X  n  matrix  function  of  t«  x  which/  for  each  x,  is  a  linear 
function  of  u(t)*  It  is  not  hard  to  show  that 

IGjt,x)l  <  Ylu(t)l/  t  e  [0,T],  X  e  C  ,  (2.41) 

where  C  is  a  compact  subset  of  0  Including  a  neighborhood  of  the  periodic 
solution  x(t)  and  y  >  positive  constant.  If  we  denote  the  fundamental 
solution  matrix  for  (2.40)  which  reduces  to  the  identity  at  t  ■  0  by 
♦^(t)/  then 

i^(t)  -  G(y<t))«^(t)  +  G^(t,y(t))*^(t) 
and  the  variation  of  parameters  formula  applies  to  give,  with  ^(t)  as  in 
(2.24), 


«^(t)  -  *(t)*^(0)  +/J  ♦(t,s)G^{s,y(s))*^(s)ds 
•  %{t)  +  ^(t,s)G^(s,y(s) ) 4^(s)ds 

from  which,  using  (2.40),  a  very  easy  argument  shows  that 

lim  4  (t)  ■  $(t)  . 

«u«  0 

Then  all  of  the  arguments  set  forth  earlier  for  (2.14)  apply  with  4(t) 
replaced  by  4  (t),  provided  lul  is  sufficiently  small,  to 

lA  1 

establish  again  the  existence  of  p  and  T  so  that  the  solution  w(t)  of 
(2.1)  with  w(0)  ■  Xq  +  p  satisfies  (2.11).  With  this  we  may  regard  the 

proof  of  Theorem  2.1  to  be  complete. 

The  next  step  is  to  establish,  rather  easily,  a  certain  localisation 
result  for  the  periodic  solutions  w(t)  resulting  from  control  functions 
u(t)  with  lul  ,  small.  Let  us  note,  as  in  (2.19)  earlier,  that, 

given  any  ^  >  0,  we  can  assume 
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ly(t)  -  x(t)  I  <  6,  t  e  [0,T^]  , 

where  y(t)  Is  now  the  solution  of  (2*38),  by  talcing  lul  -  <  C/  e 

sufficiently  small*  Then,  since  x(0)  *  x(Tq),  we  can  further  assume,  given 
&  >  0,  that 

ly(0)  -  y(T)  I<A,Tjj-t<T<Tq  +  t  ,  (2.42) 

provided  both  e  and  t  are  taken  sufficiently  small. 

How  consider  the  map  Z(p,T)  «  s(p,T}  -  p  with  z(p,t)  ■  z(t)  the 
solution  of  (2.39)  satisfying 

z(0)  •  p  e  P  . 

Since  3(z(p,T)  -  p)/3p  Is  the  restriction  of  the  oiap  4^(T)  "1  to  P,  and 
since  p.|  e  P  Is  (modulo  scalar  multiplications)  the  only  null  vector  of 
4(T)  -  X,  It  follows  that  there  are  positive  numbers  p  and  p  such  that 
for  lul  9  <  £/ 

lz(T)  -  z(0)l  «  lz(p,T)  -  pi  >  yipl 
for  Ip  I  <  2p,  Tq  -  T  <  T  <  Tg  +  T.  Then,  by  taking 

A  <  pp 

In  (2.42)  (which  may  Involve  further  restriction  of  e  and  t)  we  see  that 
we  cannot  have 

y(T)  +  z(T)  ■  w(T)  ■  w(0)  •  y(0)  +  z(0)  (2.43) 

for  z(0)  "Xg  +  p,  pep,  p<  Ipl  <  2p,  since  (2.43)  would  give 

y(T)  -  y(0)  -  z(0)  -  z(T) 

which  Is  Impossible  since 

ly(T)  -  y(0)l  <  A  <  UP  <  plpl  <  lz(T)  >  z(0)l 

Cor  such  p. 
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Equally  %fell,  the  proof  of  Theorem  2.1  shorn  that  for  small  Ipl/  p  6  P, 

and  appropriately  selected  T  >  0,  no  solution  w(t)  originating  at 

w(0)  ■  Xg  +  p  and  corresponding  to  a  control  u  with  lul  <  e,  e 

n  1 

sufficiently  small,  can  have  a  period  T  with  T  in  the  ranges  T^  -  2^  <  T 
<  Tg  -  T,  Tq  +  T  <  T  <  Tq  +  2T,  We  can  put  all  of  this  together  in  the 
following  "localization"  result. 


Proposition  2.2.  If  e  is  sufficiently  small  we  can  find  positive  t  and 


P  such  that  when  lul 


m  I 


(i)  equation  (2.1)  has  a  solution  w(t)  periodic  with  period  T  and 
with  w(0)  »  x(0)  +  p, 

{p,T)  e  N(P,T)  -  {pep,  Ipl,  Tjj  -  T  <  T  <  Tjj  +  t} 

(ii)  equation  (2.1)  has  no  solutions  of  period  T  and  initial  state 

W(0)  “  Xg  +  p  1^ 

(p,T)  e  H(2P,2T)  -  ll(P,T)  . 

The  localization  result  is  quite  important  in  the  study  of  the  optimal 

control  problem  Introduced  below  because  it  enables  us  to  single  out  a  set  of 

triples  X,  u,  T  which  may  be  confined  to  a  bounded  region  simply  by  talcing 

lul  -  sufficiently  small,  rather  than  by  introducing  further, 

L^[0,T] 

extraneous,  constraints  into  the  problem. 


It  seems  fairly  clear  that  Theorem  2.1  and  Proposition  2.2  can  be 

3f 

extended  to  the  general  system  (1.11)  provided  that  (x,u)  remains 
uniformly  bounded  for  x  6  C  C  0,  u  e  E™,  where  C  is  a  compact  subset  of 
e"  containing  the  trajectory  x(t),  0  <  t  <  T^j,  in  its  interior. 
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W*  proe«ed/  now,  to  Introdueo  and  study  the  optimal  control  problem 
referred  to  in  Section  1.  For  the  mooent  we  replace  the  general  cost 
Integrand  w(x,u)  of  (1.12)  by  a  less  general  expression/  augmented  by  a 
parameter  o, 

* 

W(  O/X)  -t-  U  OU  * 

Zt  should  not  be  difficult  to  extend  the  work  to  more  general  integrands 
w(  C/X/u) . 

Ziet  the  system  (2.1)/  i.e., 

X  -  g(x)  +  H(x)U/  X  6  s'*/  u  6  B® 

have  the  properties  set  forth  preparatory  to  and  in  Theorem  2«1  and  let  x(t) 
be  a  periodic  solution  with  least  positive  period  Tq  of 

X  -  g(x) 

% 

having  the  properties  developed  earlier.  I.et  U  be  a  positive  definite 
m  X  m  matrix  and  let  W(  O/X)  be  a  continuous  function  of  O/X  for  o  >  0 
and  X  e  0  C  r”  such  that 


W(0,x(t))  = 

0/  t  e  [O/Tq]  , 

(2.44) 

W(  0/X)  >  0, 

0  >  0/  X  S  0  / 

(2.45) 

W(  0/X(t)  ) 

4  0/  o>  0/  t  e  [0,Tq]  . 

(2.46) 

For  each  trajectory,  control  and  period  triple  w,  u,  T  with  w(t),  u(t) 
defined  on  the  interval  [0,T.j]/  and  w(0)  “  w(T),  0  <  T  <  T^,  we  define  the 
cost  functional 

J(0/W,u,T)  “  IJ  CW(0/W(t))  +  u(t)  OU(t)]dt  .  (2.47) 

We  require  >  Tq  and  we  denote  by  M  the  sot  of  trajectory,  control  and 

2 

period  triples  w.  u.  T.  defined  for  0  <  t  <  T,,  u  6  L_[0tT.] .  such  that 

I  n  1 

w,  u  satisfy  (2.1)  on  [O.T^] .  0  <  T  and 

w(0)  ■  w(T)  . 
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There  may  be  several,  or  even  infinitely  many  values  of  T  corresponding  to  a 
given  pair  w,  u.  Theorem  2.1  shows  that,  given  the  existence  of  the  periodic 
solution  x(t),  with  period  T»,  corresponding  to  the  control  u(t)  =  0 


i 


with  appropriate  assumptions  on  ♦(T^),  M  includes  at  least  one  triple  w,  u. 


T  for  every  u  6  L  tO,T^]  with  <ul  _ 

in  1 


sufficiently  small.  We  make  M 


into  a  metric  space  by  defining  the  distance  function 

^  ^  ^  ^  ^ 
d(w,u,T;  w,u,T)  »  d  (R(w),R(w))  +  lu-ul 

”  L^0,T^] 


+  T  -  T 


where  R(w)  *  (w  =  I  w  *  w(t)  for  some  t  6  [0,T^]}  and  d^  is  the  usual 
Haussdorff  metric  on  compact  sets. 

Property  (2.46)  of  the  function  W  shows  that  if  o  >  0 

J(o,x,0,Tg)  2  Jjj(O)  >  0  , 

and  the  continuity  of  W  with  respect  to  o  shows  that 

11m  J  (O)  *  0  . 

a*o  ” 

For  o  >  0  we  define  to  be  the  subset  of  M  for  which 

o 

J(o,w,u,T)  <  . 

If  w,  u  T  e  M^,  it  is  clear  that 


which  Implies  that 


^  u(t)*Uu(t)dt  <  Jq(0)  , 


2  Vo'®’ 

dt  <  - -  <  — - 

1  1 


(2.48) 


(2.49) 


when  is  the  smallest  eigenvalue  of  the  symmetric  positive  definite 

matrix  U. 


The  localization  result  expressed  in  Proposition  2.2  shows  that  the  set 


( 


•J 


of  all  w,  u,  T  in  M  for  which  <  e  nay,  if  e  is 

\tO,T^] 

sufficiently  small,  be  decomposed  in  a  natural  way  into  two  components!  those 

lying  in  N(p, t)  and  those  lying  outside  N(2p,2t),  there  being  a  triple  w, 

2 

u,  T  in  N(p,t)  for  every  u  in  L  [0,T,]  with  lul  _  <  e.  The 

■  ’ 

n  1 

importance  of  this  result  lies  in  the  fact  that  it  provides  us  with  a  priori 

bounds  on  a  certain  subset  of  the  periodic  orbits  associated  with 
2 

u  e  L^[0,T^]  of  small  nora  without  the  necessity  of  imposing  any  additional 
constraints. 

Theorem ^^2. 3.  Let  e  >  0  be  sufficiently  small  amd  let  p,  t  be  such 
that  Proposition  2.2  applied  for  I'll  5  <  e*  Let  o  >  0  be  small 

l.^[0,Ti] 

enough  so  that  (cf.  (2.49)) 


Wit 

JqCo)  <  ^  .  (2.50) 

A  A  ^ 

Then  there  exists  at  least  one  element  w,  u,  T  ^  N(p, t)  such  that 

J(a,w,u,T)  <  J(a,w,u,T),  w,u,  T  6  N(p,T)  . 

Proof .  The  Caratheodory  regularity  theory  referred  to  earlier  allows  us  to 

see  that  if  e  >  0  is  sufficiently  small  and  w,  u,  T  is  in  N(p, x)  with 

Hu  I  _  <  c,  then  there  exists  B( e)  >0  such  that 

IQ  1 


Mt)  -  x(t)  I  <  B(e),  t  e  [0,T^]  . 

Moreover , 

lim  B{  e)  “  0  . 

c-K) 

In  particular,  then,  there  is  some  M  >  0  and  a  closed  subset 
C  C  O  n  {w  I  llwJ  <  M}  such  that 

w(t)  e  c,  t  e  [o,T^]  , 


4 
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for  such  w,  u,  T.  Using  the  Schwarz  Inequality  we  quickly  see  that 


where 


Mt)  -  w(t+5t)*  -  (g(w{s))  +  H(w(s))u(a))ds* 

<  g^5t  +  h^  lu(a)*da 

<  t  h,  •u(s)|2ds)'/2  (r*"  1  ds)‘ 

I  1  u  ^ 


<  cf-^t  +  6t 


g  *  sup  Ig  (w) <  , 

'  wee 


h  ■  sup  *H(w) <  • 

'  wee 


since  (2.50)  is  true,  (2.49)  shows  that  we  nay  assume  a  minimizing  sequence 
Ujj,  Tj^  such  that  (starting  with  k  =  2  to  avoid  confusion  with 


already  defined) 


Wk,  u^,  e  N(P,T),  k  -  2,3,***, 


T.J„(0) 


2 

tn  1 


<  e,  k  -  2,3,*»*< 


lim  J(cr,w  ,u  ,T  )  -  inf 
k-*-»  K  K  K  w,u,TeN(P,T) 


inf  J(a,w,u,T)  .  (2.51) 


Prom  the  boundedness  and  equi continuity  of  the  w^,  the  boundedness  of  the 

2  ...  . 

T,.  and  the  weak  compactness  of  (u  6  I."tO,T  ]  1  lul  _  <  — ~  )  we 

m  1 

conclude  that  we  can  find  w,  u,  T  in  N(P,T)  and  a  subsequence  of 
{w^,u^,T^)  which,  for  convenience,  we  will  still  call  {wj^,u^,T^}  such  that 


# 


(2.52) 


lim  (sup  Iw.(t)  -  w(t)l]  -  0 
k+«  teT^ 


llB  T  -  T  (2.53) 

k-*« 

and  {Uj^}  converges  weakly  to  u  in  L^{0,T] .  Combining  (2.52)  and  (2.53) 
it  is  clear  that 


T  A 

lln  ^  /g’'  W(o,W|^(t))dt  -  ^  i?  W(o,w(t))dt  . 

k^*  k  T 


From  the  weak  convergence  of  {«l  }  to  u  in  L^(0,T.]  we  immediately  have 

ml 

the  weak  convergence  of  {u^^}  to  u  in  I>^I0,T).  Then,  since  for  any 

V  e  L^IO.T,] 


T  T 

I/a’^  V(t)\(t)dt!  <  (4’'  lu^(t)l^dt)^^2 

(/I^  «v(t)|2dt)’/^  <  elvl 


(or  L^[T.  ,T]  if 

tn  K  . 


<  T)  and,  as  is  well  known 


11m 

k 


f(t)dt  *  0 


for  any  Integrable  f(t),  we  conclude  that  if  the  u^(t)  are  redefined  to  lae 
zero  in  tTjj,T^],  k  *  2,3,  and  u(t)  is  redefined  to  be  zero  in  [T,T^], 
it  will  stll,l  be  true  that  (u^)  converges  weakly  to  cl.  Then,  using  a 
result  in  [H]  (Theorem  8,  p.  209)  we  know  that 


^  *  T 

/J  u(t)  Uu(t)dt  <  lim  Uj^(t)*Uu^(t)dt 
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and  we  see  that 


J(  a»w,u,T) 


<  lim  J(o/Wj^,Uj^,T  ) 
k-Hi* 


and  hence,  from  (2.51),  that 

A  A  A 

J(o,w,u,T)  ••  inf  J(o,w,u,T) 

w,u,TeN(  p,  t) 


«  lim  J{  o,W|^,Uj^,Tj^) 
k+«» 


and  the  theorem  is  proved. 

It  will  be  seen  in  the  course  of  the  discussion  in  Section  3  that  the 
above  optimal  control  formulation  is  not  adequate  for  all  purposes.  In  some 
cases  we  need  to  take  a  fixed  cost  functional 

J(x,u,T)  -  ^  lW{x(t))  +  u(t)*Uu(t)]dt 

where,  assuming  the  system  (2.1)  now  augmented,  as  in  Section  1,  with  a 
parameter  6,  viz.: 

x-f(x, e,u)  (2.54) 

or 

X  *  g(x,9)  +  H(x,  9)u  ,  (2.55) 

we  suppose  that  W(x(t))  50  on  a  T-periodic  solution  x(t)  of  the 
"uncontrolled"  system  (2.54)  or  (2.55)  with  u  =  0  and  9*9,  where  9  is 
some  minimal  value  of  the  parameter  p,  and  that 

W(x)  >  0,  X  ^  {x(t)  I  0  <  t  <  T}  . 

The  optimal  solution  for  9*9  is  clearly  x(t),  0  <  t  <  T.  One  wishes  to 
modify  Theorem  2.3  to  establish  the  existence,  locally,  of  an  optimal  solution 
for  9  in  some  neighborhood  of  9. 
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The  only  modifications  that  are  at  all  significant  here  involve 
modification  of  the  localization  result.  Proposition  2.2,  so  that  it  applies 
for  lul  -  sufficiently  small  and  6  sufficiently  near  6  and 

m  1 

replacement  of  the  ^  priori  bound  ( 2 . 48 )  on  u  by 

T 

i  fl  u(t)*Uu(t)dt  /q*'  W{x  (t))dt 

w 

where  x  (t)  is  some  periodic  solution  of 

y 

X  (t)  -  f(x  (t),p,0} 

y  y 

lying  near  the  periodic  solution  x(t)  with  period  near  T  -  the 

existence  of  which  is  established  much  as  in  Theorem  2.1,  except  that  the 

proof  can  be  greatly  simplified  now,  using  the  implicit  function  theorem 

since,  regarded  as  a  control,  y  is  now  (trivially)  continuous  in  t.  In 

fact,  X  and  T  are  unique  in  this  case.  The  rest  of  the  argument  is 

y  y 

almost  word  for  word  the  same  as  in  Theorem  2.3  above  2md  will  not  be  repeated 
here. 
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3.  N»c»««ary  Conditions  for 


ZZSnSBSTU 


Coaputatlonal  Conaidarationg 


Xn  the  preaent  section  %ie  consider  the  system 

X  -  f(x,u),  X  e  O  c,  b",  u  e  E®  ,  (3.1) 

wherein  It  Is  assumed  that  f(x,u)  Is  continuously  differentiable  with 
respect  to  both  x  and  u  over  the  range  Indicated.  We  suppose  that  x(t) 

A 

Is  a  periodic  solution  of  (3.1),  with  least  positive  period  T,  which 

A  ^ 

corresponds  to  the  control  u(t)  on  [0,T1.  For  the  moment  we  malce  no 
further  assumptions  on  f  but  restrict  attention  to  piecewise  continuous  u 
and  piecewise  continuous  variations  ^u  from  u,  but  fairly  minor 
modifications  allow  us,  in  the  case  of 

X  ■  g(x)  +  H(x)u,  X  e  0  C  e”,  u  6  E®  ,  (3.2) 

the  system  treated  In  Section  2,  to  also  admit  the  possibility  that  u  and 
the  variations  ^u  are  just  square  Integrable. 

h  number  of  authors  have  considered  necessary  and  sufficient  conditions 
for  periodic  systems  ([I],  [J] ,  [K] ,  (L]).  We  do  not  develop  the  conditions 
here  t>ecause  they  are  new  -  they  are  not  **  but  to  allow  us  to  make  appropriate 
reference  later  in  this  section.  We  seek  necessary  conditions  in  order  that 
the  control  u  should  afford  a  local  minimum  (at  least)  for  the  cost 
functional 

J(x,uT)  «  v(x(0))  *  ^  w(x(t) ,u(t) )dt  ,  (3.3) 

where  v(x)  and  w(x,u)  are  continuously  differentiable  in  x,  (x,u) 
respectively,  for  x  e  O  c  e",  u  6  b"'.  In  (3.3)  it  is  assumed  that  x(t)  is 

A 

a  periodic  solution  of  (3.1)  having  period  near  the  period  T  for  x(t)t 

It  -  ^1  <  T,  T  >  0  , 

that  x(t)  lies  near  x(t)  in  the  sense  that 

»x(t)  -  x(t)»  <  e,  0  <  t  <  sup(T,T) 
for  some  e  >  0  and  that 
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Iu(t)  -  u(t)l  <  6,  0  <  t  <  aupCT^T)  , 

for  aoM  6  >  0*  In  tha  eaaa  of  (3>2)  this  can  ba  ralaxad  to 

lu  -  ul  .  <  5/  T.  >  ■ax{T,T)  • 

L^10,T  ]  ^ 

Lat  6a(t)  ba  piaoawlaa  oontlnoous  for  0  ^  ^  conaldar 

oontrola 

a(t)  -  tt(t)  a  fla(t),  a  raal  ,  (3*4) 

^  A 

whara  a(t)  la  axtandad  outalda  [0#T]  by  parlodiclty  as  raqulrad* 

Wa  su^osa  tte  tlpa  paraa«tar  normaliaad  in  suoh  a  way  that  u(t)  is 

A 

continuous  at  t  >  0  (aquivalantly  at  T)  and  wa  sat 

F(t)  -  F(x(t),  u(t))  -  ~  {x(t),  u(t)>»  0  <  t  <  T  •  (3.5) 

hs  aarliar  in  tha  papar,  wa  danote  by  4(t,s)  tha  Mtrix  solution  of 

-  F(t)*(t,o),  «(s.s)  -  1  , 

and  wa  abbraviata  4(t«0)  by  4(t). 

Tha  variational  thaory  for  a  controllad  pariodic  solution  x(t)  is 
sooMwhat  difforant  from  that  for  an  uncontrollad  pariodic  solution  x(t)  as 
discussad  aarliar  in  this  papar.  If  wa  lat 

p  ■  f(x(0),  a(0))  (3.6) 

A 

it  is  no  longar  nacassarily  trua  that  p  it  an  aiganvaotor  of  4(T) 
corrasponding  to  an  aiganvalua  ■  1.  To  saa  why  this  is  so.  contidar. 
a.g./  tha  casa  whara  u(t)  is  pariodic  and  of  class  on  [O.T].  Lat 

p(t)  ■  x(t)  •  f(x(t).  u(t)).  0  <  t  <  T  . 

Oiffarantiating  (3.1)  givas 

p  -  F(t)p  ♦  H(t)u(t)  ,  (3.7) 

whara 

H(t)  «  H(x(t),  u(t))  ■  ^  (x(t),  U(t))  ,  (3.8) 


so  that 


p(t)  -  ♦(t)p(O)  +/J  «(t,s)H(s)tt(s)ds  . 

A  A  A  A 

Since  p(T)  <■  x(T)  -  x(0)  >  p(0)  «  p,  we  have 

A  • 

(I  -  ♦(T))p-  /J  ♦(T,8>H(8>u(8)da  (3.9) 

A 

and  p  is  an  eigenvector  of  «(T)  corresponding  to  the  eigenvalue  1  If  and 
only  If  the  right  hand  side  of  (3*9)  turns  out  to  be  sero*  Of  course^  In  the 
uncontrolled  case  the  right  hand  side  Is  automatically  equal  to  zero. 

He  are  going  to  proceed  under  the  following 
Assumption >  Either  (Dt  p,  given  by  (3*6),  Is  the  unique  eigenvector  of 

A  A 

i(T)  corresponding  to  the  eigenvalue  1,  which  Is  simple,  or  (ll)t  6(T) 
has  no  eigenvalue  equal  to  1 . 

A 

Zt  Is  not  particularly  difficult  to  treat  the  case  where  6(T)  may  have 
a  simple  eigenvalue  1  with  corresponding  eigenvector  different  from  p  but 
we  shall  not  do  so  here. 

In  case  (1)  the  study  of  perturbed  periodic  solutions  corresponding  to  a 
small  perturbation  5u(t)  from  the  control  u(t)  Is  just  Theorem  2.1.  In 
case  (11)  one  can  establish  the  following 

A 

Prqgos it Iot  3 . 1 »  If  ♦(T)  has  no  eigenvalue  equal  to  1  then, 
corresponding  to  each  piecewise  continuous  perturbed  control 

u(t)  “  u(t)  +  o  6u(t) 

with  I  a|  sufficiently  small,  there  is  a  unique  perturbation  6x(0)  of  the 
Initial  state  so  that  the  solution  x(t)  o^  (3.1)  corresponding  to  the 
control  u(t)  and  the  initial  state  x(0)  -  x{0)  +  Sx(0)  remains  periodic 

A 

with  period  T. 

He  omit  the  proof.  A  comparable  result  can  be  obtained  for  (3.2)  with 
5u  in  L^I0,T]. 
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To  procoed,  we  consider  a  perturbed  control  (3*4),  %fhich  «fe  take  to  be 


piecewise  continuous  in  this  discussion,  and  we  let 

x(t,a)  ■  x(t)  +  a  ix(t)  +  o(|oi|),  |o|  ♦  o 


(3.10) 


be  a  perturbed  periodic  solution  near  x(t)  (uniqueness  falls  under  ease  (ii) 
of  our  Assumption  even  with  piecewise  continuous  controls)  with  period 

T(a)  -  T  +  o  fir  +  o(|aj),  |o|  ♦  0  . 

In  (3.10)  Sx(t)  is  a  solution  of  the  variational  equation  (cf.  (3.5),  (3.7), 


(3.8)) 


fix(t)  -  P(t)6x(t)  +  H(t)fiu(t) 


(3.11) 


for  0  <  t  <  T,  and  the  equation  la  extended  outside  this  range,  as  required, 

using  the  T  periodicity  of  x(t)  and  u(t). 

Denoting  the  gradients  of  v,  w  with  respect  to  x,  (x,u), 

9v  3w  3w 

respectively,  by  g^,  and  treating  these  as  row  vectors,  it  may  be 

seen  that 

J(x(»,o),  u(*,a),  T(a))  -  j(x,u,T) 


+  {|“  (?i(0))«x(0)  +  ^  /f  l|5  (x(t),  u(t))6x(t) 

T 

3  A 

+  g^  (x(t),  u(t))6u(t)]dt  w(x(t),  u(t))dt 


+  —  w(x(0) ,  u(0) ) )6t}  a+o(|a|),  o*0 


(3.12) 


wherein  we  have  used  periodicity  to  see  that 

w(x(T),  u(T))  -  w(x(0),  u(0))  . 

Now  the  variational  derivatives  6x  and  are  not  independent* 

the  last  two  are  determined  by  ^u  and  the  normalization  of  5x(0).  In  order 
to  be  able  to  express  the  cost  variational  derivative  explicitly  in  terms  of 


32- 


■  ■  X.  -  ww; 


Au,  rather  than  Implicitly  through  5x  and  one  introducea  a  solution 

l(t)  of  the  "adjoint  system" 


X(t)  -  -P(t)*X(t)  +  ^  (x(t),  u(t))  .  ( 

T 

The  boundary  conditions  will  be  specified  shortly*  We  find  then  that 


X(T)*«x(T)  -  X(0)*6x(0)  -  /J  ~(X(t)*6x{t))dt 


A 

(of.  (3.11),  (3.13))  -  /T  I(-X(t)*P(t)  U(t))«x(t) 


+  X(t)  (P(t)«x(t)  +  H(t)«u(t))]dt 


(3.13) 


A  ^  A  . 

X(t)  H(t)6u(t)dt  A  3x  u(t))fix(t)dt  . 

T 


To  first  order,  the  periodicity  requirement 


(3.14) 


x(T  +  5t)  ■>3^(0) 


may  be  replaced  by  the  condition 


and  then,  since 


A  A 

5x(T)  +  -r-  (T)6t  ■  fix(O) 
at 


~  (T)  ■  f(x(T),  u(T)) 


we  have 


f(x(0),  u(0)) 


«x(T)  +  f(x(0),  u(0))fiT  -  «x(0)  . 


(3.15) 


Using  (3*14)  and  (3*15)  in  (3*12)  we  see  that 


•xW”. 


^  (xC,ci),  uC.o),  T(o))l^^  - 

-  ^(0)*)«x(0)  -  {X(T)*f(x(0),  u(0)) 

A  " 

-  —  w(x(0),  u(0))  +  ^  %»(x(t),  u(t))dt)iT 

T 

+  /:  txU  (x(t),  u(t))  -  X(t)  H(t)]«u(t)dt  .  (3.16) 

The  dependence  bii  fix(O)  le  annulled  by  imposing  the  boundary  condition 

X(T)  -  X(0)  +  -g  (x(0))  -0  .  (3.17) 

At  this  point  vfe  have  to  consider  two  possibilities,  corresponding  to  the 

A 

t%iK>  cases  delineated  in  the  above  Assumption.  If  case  (i)  obtains,  *(T)  has 
a  simple  eigenvalue  of  1.  Then  the  fundamental  solution  matrix 

’f(t)  -  (♦(t)*)"^  , 

obeying 

-F(t)*'r(t) 

tM  ^ 

l8  such  that  t(t)  has  an  eigenvalue  equal  to  1  and 

X(t)  -  -F(t)*X(t) 

has  a  one  parameter  family  of  T-periodic  solutions  which  tre  may  represent  in 

A  . 

the  form  0X(t).  Given  a  solution  X(t)  of  (3. 13), . (3. 17) 

X(t)  -  X(t)  +  eX(t) 

also  satisfies  (3.13),  (3.17)  for  every  value  of  0.  Then  tfe  may  expect  to  be 
able  to  select  0  (with  suitable  non- degeneracy  requirements)  so  that  X(t) 
also  satisfies 
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X(T)  f(x(0),  u(0))  -  ^  w(x(0),  u(0)  *  it  w(*(t),  u(t))dt  *0  .  (3.18) 

T 


&2  '0 
T 


With  A(t)  so  chosen,  we  see  from  (3.16)  that  d7/da  >  0  for  all  ju(t) 
just  in  case 

X(t)*H(t)  -  ^1“  (x(t),  u(t))  -  0,  t  e  [0,T]  .  (3.19) 

T 

If  case  (11)  obtains  we  have  noted  that  we  can  find  perturbed  periodic 
solutions  for  which  5t  *  0.  This  also  annuls  the  second  term  In  (3.16)  and 
enables  us  to  conclude  that  (3.19)  must  hold.  But,  In  ease  (11)  still.  It  Is 
also  possible,  for  each  5u  and  a  sufficiently  small  in  (3.4),  to  find 
perturbed  periodic  solutions  x(t,a)  near  x(t)  with  periodic  T  -f  a  6t, 

0  given  In  advance  (this  is  true  because  (cf.  (3.15))  we  can  always 

solve 

(♦(T)  -  l)«x(0)  +  f(x(0),  u(0))«T  -  0 

In  ease  (ID).  Assuming  (3.19)  already  established  as  noted,  dj/do  reduces 
to  the  second  term  in  (3.16)  and  we  conclude  that  condition  (3.18)  must.  In 
fact,  hold  in  this  case  also. 

It  will  be  observed  that  if  (x(0))  *  0,  (3.17)  becomes 

ox 

X(T)  -  X(0)  -  0 

,  A 

SO  that  A(t)  is  periodic  with  period  T. 

It  is  customary  to  refer  to 


H{X(t),  x(t),  u(t),  T) 

-  X(t)*f(x(t),  u(t))  -^w(x(t),  u(t)) 


(3.20) 


as  the  Hamiltonian  for  this  optimization  problem.  The  condition  (3.19)  is 


2h  A,.  .  A,  .  A.  .  ..A. 


3u 


(X(t),  x(t),  u(t),  T)  -  0,  t  e  [0,T]  . 


(3.21) 


We  have  the  following  result. 
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Proposition  3»1:  H(X(t)«  x(t),  u(t),  T)  Is  constant  for  t  6  [0,T] 

Proof.  Pron  (3.5),  (3.8),  (3.20)  we  have 


~H(X(t),  x(t),  u(t),  T)  -  X(t)*£(x(t),  u(t)) 


+  (X(t)  F(t)  “^•^(x(t),  u(t)))x(t) 
T 


+  (X(t)  H(t)  "  ^  U(t)))u(t) 

T 


(using  (3.19),  (3.1),  (3.13)) 


(-X(t)*F(t)  +  X  (x(t),  u(t))]x(t) 
T 


+  [X(t)  F(t)  “  X  (x(t),  u(t))Ix(t)  »  0  . 

T 

When  X(T)  -  X(0),  (3.19)  will  ordinarily  give  tl(T)  -  u(0)  and, 

x(T)  *  3^(0),  (3.18)  becomes 

A 

H(X(t),  3(t),  u(t),  ^  w(x(a),  u(s))ds,  t  e  [0,T]  .  (3. 


If  we  let 


U(t) 


^  X(t), 


t  8  10,T)  , 


we  have,  for  t  8  I0,T] , 


C(t)  -  -F(t)*M(t)  +  -^  (x(t),  u(t))  , 


(3, 


U(t)*H(t)  -  ~  {x(t),  u(t))  =  0 


(3. 


where 


•^H(M(t),  x(t),  ti(t))  -  0  , 


(3. 


H(p(t),  x(t),  u(t))  *  y(t)  f(x(t),  u(t))  -  w(x(t),  u(t))  .  (3. 


since 

21) 

22) 

23) 

24) 

25) 
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(3.26) 


In  this  case 

J(x,u,T)  =  H(U,(t),x(t),u(t))  +  v(x{0)),  t  e  [0,T]  . 

Very  little  needs  to  be  changed  in  the  development  which  leads  up  to 
Theorem  3.2  in  order  to  permit  u  to  be  a  function  L^CO.T]  in  the  case  of 
the  special  class  of  systems  (3.2).  In  fact,  the  class  is  selected  so  that 
with  controls  u(t,a)  as  in  (3.4),  but  u(t),  6u(t)  now  in  L^[0,Tl,  the 
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resulting  solution  x(t,a)  still  has  the  form  (3«10),  with  5x(t)  satisfying 

(3<11).  For  general  functions  f(x,u)  It  Is  not  easy  to  see  that  this  Is  the 

case  without  making  some  rather  technical  assumptions  on  f  to  account  for 

the  possibility  that,  while  I6ul  in  L^[0,T]  might  be  small,  5u(t)  might 

still  assume  arbitrarily  large  values.  The  form  (3.2)  can  be  replaced  by 

3f 

various  boundedness  assumptions  on  but  we  elect  not  to  pursue  this  here. 

A  2  ^ 

The  time  variable  In  the  case  of  (3.2)  with  u  In  L  [0,T]  should  be 
normalized  In  such  a  way  that  0  Is  a  Lebesgue  point  for  u,  with  the  effect 

-ft 

that  (0)  Is,  indeed,  equal  to  9(x(0)  ■<-  H(x(0))il(0),  so  that  (3.15) 
appropriately  modified,  still  makes  sense.  (Of  course  u(t)  turns  out  in 
nearly  all  cases  to  be  smooth  and  periodic.) 

It  is  not  generally  very  easy  to  prove  the  existence  of  solutions  of  the 
equations  constituting  this  set  of  necessary  conditions.  However,  under  the 
circumstances  %rtiich  correspond  to  Theorem  2.3  it  is  possible  to  p^ove  the 
existence  of  a  unique  solution  (In  a  certain  sense)  for  small  values  of  the 
parauneter  0,  provided  we  adjoin  to  (2.44),  (2.45),  (2.46)  the  further 
condition  that 

(0,x(t))  =0,  t  e  [0,Tjj]  (3.30) 

and  the  function  v(x)  appearing  in  (3.3)  satisfies  certain  requirements, 
about  which  we  will  have  more  to  say  below.  Let  us  notice  first  of  all  that 
when  v(x)  =  0  there  is  no  possibility  of  a  unique  minimum  for  the  cost 
functional  (3.3)  under  the  stipulated  circumstances  l>ecause  translation  by  a 
fixed  5t  along  the  periodic  optimal  trajectory  x(t),  li(t)  ((3.27)  shows 
M(t)  to  be  periodic  when  v(x)  =0)  yields  a  new  trajectory  with  the  same 
cost  as  well  as  a  new  set,  strictly  speaking,  of  solutions  of  the  necessary 
conditions.  This  sort  of  non-uniqueness  Is  not  very  interesting,  of  course. 
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and  It  Is  customary  and  convenient  to  rule  It  out  by  Imposing  an  appropriate 
normalizing  condition  on  the  initial  state  x(0).  Letting  (cf*  Section  1) 

p  -  x(0)  -  f(x(0),0) 

one  can  require,  for  example,  that  a  nearby  optimal  x(t)  should  satisfy 


p*(5(0)  -  x(0))  -  0 


(3.31) 


thereby  confining  the  initial  state  ^(0)  to  an  (n-’l ) -dimensional  hyperplane 
transverse  to  the  nominal  orbit  x(t).  Computationally  it  is  preferable  not 


to  enforce  (3.31)  as  an  explicit  constraint  but,  instead,  to  use  for  v(x)  a 


function  such  as 


v{x)  -  (p*(x(0)  -  x(0)))^ 

=-  (x(0)  -  x(0))*pp*(x(0)  -  x(0))  . 


(3.32) 


In  obtaining  the  least  possible  cost  for  (3.3)  one  automatically  obtains  the 
periodic  x(t)  which  satisfies  (3.31)  emd,  in  the  end,  —  (x(0))  *0  so 
that  (3.27)  applies  to  show  u(t)  to  be  periodic. 


Provided  a  correctly  constructed  normalizing  v(x)  is  used  as  indicated 
above,  it  is  possible  to  use  the  implicit  function  theorem  to  construct 
x(t),  u(t),  near  x(t),  0,  satisfying  the  necessary  conditions  set  forth  in 
Theorem  3.2.  One  takes  x(0),  p{0)  to  lie  in  a  small  neighborhood  of  x(0), 
0.  Then,  via  solution  of  the  differential  equations  (3.1),  (3.13)  with  u(t) 


determined  by  (3.19),  the  vectors 


A  A  A 

x(T)  -  x(0)  , 


(3.33) 


A  A  Av  A 

W(T)  -  W(0)  +  T  (X(0)) 


(3.34) 


and  the  scalar  quantity  (cf.  (3.28)) 


H(  (j(0)  ,x(0)  ,u(0) )  +  v(x(0))  -  7  w(x(t)  ,u(t))dt  (3.35) 

T 
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I 


are  seen  to  constitute  a  (2n+1)  dimensional  vector  function  of  x(0),  M(0), 

T  and  any  relevant  parameters,  such  as  o  of  Theorem  2.3  or  the  parameter 
V,  as  in  the  discussion  of  the  alternate  problem  following  Theorem  2.3.  The 
Jacobian  of  (3.33),  (3.34),  (3  35)  with  respect  to  x(0),  )i(0)  and  T  can  be 
computed  with  the  aid  of  the  variational  equations  and  under  appropriate 
clrciimstances  (e.g.  the  cost  functional  has  the  form  discussed  in  Theorem  2.3, 
(3.30)  holds  and  v(x)  is  defined  as  in  (3.32))  one  ceai  show  that  this 
Jacobian  is  nonsingular  when  0^0  (or,  depending  on  the  problem,  when 
P  =  0 )  and  unique  x(0),  P(0),  T  can  then  be  Inferred  for  o  (or  V) 
sufficiently  small.  We  will  not  go  into  detail  on  this  as  it  will  appear 
elsewhere  and  would  lead  us  much  too  far  astray  here. 

The  computational  version  of  the  above  is  the  simple  Newton's  method  for 
solving  the  equations  obtained  by  setting  (3.33),  (3.34),  (3.35)  equal  to 
zero.  Implementation  requires  numerical  integration  of  the  variational 
equations  over  an  interval  10,T^],  where  T^  is  the  latest  approximation  to 
the  period  ^  in  order  to  be  able  to  compute  the  Jacobian  of  (3.33),  (3.34), 
(3.35)  with  respect  to  the  current  initial  conditions  for  Xj^,  and  the 
period.  Just  as  noted  above,  it  can  be  established  that  this  procedure  is 
locally  convergent  in  certain  cases.  However,  this  method,  while 
theoretically  useful  and  appealing  in  its  simplicity,  suffers  from  a  number  of 
disadvantages  as  a  computational  procedure.  The  necessity  of  integrating  the 
variational  equations  to  obtain  the  Jacobian  matrix  is  one  such  disadvantage 
to  begin  with,  but  more  serious  is  the  frequent  numerical  instability  of  the 
procedure.  In  most  cases  the  periodic  solution  x(t),  corresponding  to 
u(t)  •  0,  that  one  begins  with,  is  asymptotically  stable.  The  corresponding 
variational  system  is  a  time  varying  periodic  system  whose  characteristic 
exponents  have  negative  real  parts.  But  the  dual  variational  system  for  the 


adjoint  variable  P  (or  X)  has  characteristic  exponents  which  are  the 
negatives  of  those  of  the  variational  system  for  x(t)*  Thus  the  variational 
system  for  V  (or  X)  has  solutions  which  may  grow  very  rapidly.  This  means 
that  It  Is  often  very  difficult  to  accurately  compute  the  part  of  the  Jacobian 
corresponding  to  the  equation  (3.34).  This  has  been  a  serious  limiting  factor 
In  our  computational  experience  thus  far. 

While  theoretically  less  rapidly  convergent,  the  steepest  descent  method 
offers  much  In  the  way  of  Improved  performance  from  the  standpoint  of 
numerical  stability,  primarily  because  the  state  equations  and  adjoint 
equations  are  Integrated  In  opposite  directions  to  take  advantage  of  their 
differing  stability  characteristics.  We  will  describe  briefly,  one  step  of 
such  a  procedure. 

We  shall  suppose  that  Xj^(t),  u^(t)  form  a  periodic  trajectory  •  control 
pair  for  (3.1)  with  period  Tj^.  As  we  do  not  require,  ab  Initio,  any 
continuity  for  Uj^,  the  periodicity  refers  to  the  trajectory  Xj^(t).  With 
x(t)  >•  Xj^(t)  available,  one  proceeds  to  compute  a  solution  (3.13). 

Corresponding  to  the  two  cases  admitted  under  the  Assumption  above,  one 
follows  two  routes.  If  case  (1)  applies,  determined  so  that  It  also 

satisfies  (3.28).  Then  setting 

5uj^(t)  -  H(t)  Mj^lt)  -  (x(t),u(t))  (3.36) 

gives  (cf.  (3.16))  dJ/da  <  0  (If  (3.23)  does  not  already  hold  for  x^,  W^). 
For  small  positive  a  we  let 

Uj^^^(t)  -  Uj^(t)  +  a  «u^(t),  t  e  [0,Tj^]  . 

Extension  of  {U|^(t)  outside  [0,T]^],  where  required,  can  be  accomplished 
from  (3.36)  and  may  assume  this  has  already  been  done  for  u^(t).  Then  we 

determine  5x^(0)  and  ^T^  so  that  the  solution  of  (3.1)  with 


initial  state  +  fiXj^(O)  and  u(t)  =  u^^^{t)  satisfies  + 

“  ’‘k+l^^k+1^  "  «  >  0  is  sufficiently  small  one  can  show 

that  the  cost  associated  with  will  be  smaller  then  that 

associated  with  x^^*  u^«  If  case  (ii)  applies  is  onl^  required  to 

satisfy  (3«27),  6U|^  is  selected  as  in  (3.36),  and  we  let  >  o  times  the 

quantity  on  the  left  hand  side  of  (3.28)i  we  will  again  see  from  (3.16)  that 
d3/da  <0.  This  defines  the  new  period  t^^^  *  ^  ^^k  ^  proceed,  as 

we  may  in  case  (ii),  to  find  corresponding  to  the  control 

having  precisely  the  period  Again  the  cost  associated  with  U|^^^ 

will  be  reduced,  ae  compared  with  Xj^,  Uj^,  if  a  >  o  is  sufficiently  small. 

The  above  procedure  is  continued  until  the  necessary  conditions  are 
satisfied.  The  integrations  necessary  to  determine  Uj^(t)  are  always  carried 
out  in  reverse  time  (as  compared  with  those  for  Xj^(t))  for  the  sake  of 
stability.  The  above  description  is  very  rou^  and  a  number  of  refinements 
remain  to  be  delineated;  for  example,  how  do  we  treat  the  situation  which 
obtains  when  case  (i)  is  very  nearly,  but  not  quite,  true.  This  is  not  the 
place  to  treat  such  questions. 

While  case  (ii)  may  be  considered  generic  for  our  problem,  case  (i)  does 
arise  in  significant  instances;  e.g.  with  a  cost  functional  (2.47),  satisfying 
(3.30),  the  initial  step  away  from  the  uncontrolled  solution  must  be  made 
under  case  (i)  circumstances. 

An  Bxample.  Initial  computational  experiments  have  been  carried  out  with  the 
second  order  scalar  system 


with  the  parameter  6  in  the  range  0  <  6  <  *6  with  the  objective  of 
amplitude  suppression.  The  uncontrolled  system  has  a  self-excited  oseillation 
with  amplitude  A(  6)  shown  by  the  solid  curve  A(  8)  in  Figure  3.2^  the  whole 
orbit  for  6  «  .6  being  sho%m  by  the  closed  solid  curve  in  Figure  3.1. 

Initially  we  worked  with  a  cost  functional  of  the  form 

^  to(x{t)^  +  (x(t))^)  +  u(t)^]dt  ,  (3.38) 

first  computing  the  uncontrolled  orbit#  which  is  optimal  for  o  "  0  and  then 
increasing  a  in  small  increments.  The  computational  method  used  was  the 
Netrton  method  described  earlier.  The  cost  functional  (3.37)  has  the  form 
discussed  in  Theorem  2.3.  It  proved  unsatisfactory#  however.  In  a  very 
narrow  range  of  positive  values  a  controlled  periodic  solution  of  smaller 
amplitude  was  obtained  but  it  "popped”#  like  a  bubble#  for  a  "  .023.  The 
word  "popped”  is  used  advisedly  because  we  believe  the  disappearance  of  t!he 
solution  near  this  value  of  a  is  indeed  due  to  mathematical  causes  similar 
to  those  which  are  operative  in  the  bursting  of  a  bubble  or  the  stretching  of 
a  soap  film}  the  stability  of  the  initial  uncontrolled  solution  manifests 
itself  initially  -  a  finite  amount  of  control  energy  is  reguired  just  to  move 
away  from  this  initial  solution.  Beyond  a  certain  point#  however#  a  minimal 
cost  solution  of  positive  amplitude  no  longer  exists  -  the  origin  is  the  only 
(degenerate)  solution.  We  are  fairly  well  satisfied  that  the  optimal  solution 
does#  in  fact,  disappear  in  "mid-air"  rather  than  collapse  continuously  to%mrd 
the  origin  -  not  only  by  numerical  results  but  also  by  related  theoretical 
considerations  which  we  shall  not  go  into  here. 

Better  success  was  achieved  by  following  the  alternate  schema  suggested 
in  the  material  which  follows  Theorem  2.3  in  Section  2.  We  used  a  fixed  cost 


functional  of  the  form 


(3.39) 


^  /J  (R(t)^  -  +  (u(t))^]dt 

where 

R{t)^  -  {x{t))^  +  (x{t))^ 

and  a  is  the  average  amplitude  of  the  uncontrolled  solution  when  9  ■  .4. 

The  condition  W(x(t))  =  0  (cf.  (2*55)  pp)  Is  not  gulte  satisfied  In  this 

case  but  the  uncontrolled  orbit  Is  so  nearly  circular  that  the  discrepancies 

have  negligible  effect.  When  6  «  .4  there  Is  no  noticeable  modifications  of 

the  orbit  but*  as  6  is  Increased  the  effect  of  the  optimal  control  policy  Is 

more  and  more  pronounced/  producing  the  modified  amplitude  curve  A  (9)  shoim 

c 

by  the  dashed  curve  In  Figure  3.2  and/  for  B  »  .6/  the  modified  orbit  shown 
as  a  dashed/  closed  curve  In  Figure  3.1. 

The  extreme  "fragility"  of  the  Mewton  method  became  apparent  In  these 
computations;  whether  the  parameter  was  o  In  (3.38)  or  9  In  (3*37)  In 
connection  with  (3*39)/  It  was  necessary  to  change  It  In  extremely  small 
Increments  to  avoid  numerical  "blow-up"  -  primarily  due  to  the  Instability  of 
the  adjoint  system  as  one  Integrates  In  the  positive  t  direction. 

In  future  work  we  hope  to  use  the  steepest  descent  method  outlined  here 
to  compute  optimal  controls  and  modified  orbits  for  a  system  representative  of 
the  wing  -  air  system  described  In  Section  1  and  shown  In  Figure  1*2.  The 
difficulties  already  noted  for  (3.38)  Indicate  that  the  use  of  the  Newton 
method  trould  be  almost  hopeless  In  this  case. 
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Figure  3.2  Controlled  and 
Uncontrolled  Amplitude  vs  Parameter 
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We  consider  a  controlled  system 

X  -  f(x,u) 

exhibiting  a  self  excited  periodic  solution  x(t)  for  u(t)  s  o  and  consider  the 
qxMstion  of  modification  of  this  orbit  to  a  controlled  periodic  solution  ft(t) 
corresponding  to  a  control  u(t)  chosen  so  as  to  minimize  a  cost  functional  of  the 
form  .  * 

I  /J  w(jl(t),  il{t))dt  , 

where  9  is  the  period  of  the  controlled  periodic  solution.  Some  relevant  applica- 
tions  are  cited. _ _ _ _ _ 
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